Scattering Theory of Gilbert Damping 
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The magnetization dynamics of a single domain ferromagnet in contact witli a thermal bath 
is studied by scattering theory. We recover the Landau-Liftshitz- Gilbert equation and express the 
effective fields and Gilbert damping tensor in terms of the scattering matrix. Dissipation of magnetic 
energy equals energy current pumped out of the system by the time-dependent magnetization, with 
separable spin-relaxation induced bulk and spin-pumping generated interface contributions. In 
linear response, our scattering theory for the Gilbert damping tensor is equivalent with the Kubo 
formalism. 
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Magnetization relaxation is a collective many-body 
phenomenon that remains intriguing despite decades of 
theoretical and experimental investigations. It is im- 
portant in topics of current interest since it determines 
the magnetization dynamics and noise in magnetic mem- 
ory devices and state-of-the-art magnetoelectronic ex- 
Sriments on current-induced magnetization dynamics 
. Magnetization relaxation is often described in terms 
of a damping torque in the phenomenological Landau- 
Lifshitz- Gilbert (LLC) equation 
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where M is the magnetization vector, 7 = g^s/fi is the 
gyromagnetic ratio in terms of the g factor and the Bohr 
magneton /is, and Ms = |M| is the saturation magneti- 
zation. Usually, the Gilbert damping (^(M) is assumed 
to be a scalar and isotropic parameter, but in general it 
is a symmetric 3x3 tensor. The LLG equation has been 
derived microscopically [3| and successfully describes the 
measured response of ferromagnetic bulk materials and 
thin films in terms of a few material-specific parameters 
that are accessible to ferromagnetic-resonance (FMR) ex- 
periments [3|. We focus in the following on small fer- 
romagnets in which the spatial degrees of freedom are 
frozen out (macrospin model). Gilbert damping pre- 
dicts a stricly linear dependence of FMR linewidts on fre- 
quency. This distinguishes it from inhomogenous broad- 
ening associated with dephasing of the global precession, 
which typically induces a weaker frequency dependence 
as well as a zero-frequency contribution. 

The effective magnetic field Hcfr — — 9F/9M is the 
derivative of the free energy F of the magnetic system 
in an external magnetic field Hoxt, including the classi- 
cal magnetic dipolar field H^- When the ferromagnet is 
part of an open system as in Fig. [Tl — 9F/9M can be 
expressed in terms of a scattering S-matrix, quite anal- 
ogous to the interlayer exchange coupling between ferro- 
magnetic layers \S,- The scattering matrix is defined in 
the space of the transport channels that connect a scat- 
tering region (the sample) to thermodynamic (left and 



FIG. 1: Schematic picture of a ferromagnet (F) in contact 
with a thermal bath via metallic normal metal leads (N). 



right) reservoirs by electric contacts that are modeled by 
ideal leads. Scattering matrices also contain information 
to describe giant magnetoresistance, spin pumping and 
spin battery, and current-induced magnetization dynam- 
ics in layered normal-metal (N)| ferromagnet (F) systems 

In the following we demonstrate that scattering the- 
ory can be also used to compute the Gilbert damping 
tensor G(M). The energy loss rate of the scattering re- 
gion can be described in terms of the time-dependent 
S-matrix. Here, we generalize the theory of adiabatic 
quantum pumping to describe dissipation in a metallic 
ferromagnet. Our idea is to evaluate the energy pump- 
ing out of the ferromagnet and to relate this to the energy 
loss of the LLG equation. We find that the Gilbert phe- 
nomenology is valid beyond the linear response regime of 
small magnetization amplitudes. The only approxima- 
tion that is necessary to derive Eq. ([T]) including G(M) 
is the (adiabatic) assumption that the frequency co of the 
magnetization dynamics is slow compared to the relevant 
internal energy scales set by the exchange splitting A. 
The LLG phenomenology works so well because hoj -^ A 
safely holds for most ferromagnets. 

Gilbert damping in transition-metal ferromagnets is 
generally believed to stem from spin-orbit interaction in 
combination with impurity scattering that transfers mag- 
netic energy to itinerant quasiparticles [3|]. The subse- 
quent drainage of the energy out of the electronic sys- 
tem, e.g. by inelastic scattering via phonons, is believed 
to be a fast process that does not limit the overall damp- 
ing. Our key assumption is adiabaticiy, meaning that 
the precession frequency goes to zero before letting the 
sample size become large. The magnetization dynam- 
ics then heats up the entire magnetic system by a tiny 



amount that escapes via the contacts. The leakage heat 
current then equals the total dissipation rate. For suf- 
ficiently large samples, bulk heat production is insensi- 
tive to the contact details and can be identified as an 
additive contribution to the total heat current that es- 
capes via the contacts. The chemical potential is set 
by the reservoirs, which means that (in the absence of 
an intentional bias) the sample is then always very close 
to equilibrium. The S-matrix expanded to linear order 
in the magnetization dynamics and the Kubo linear re- 
sponse formalisms should give identical results, which we 
will explicitly demonstrate. The role of the infinitesi- 
mal inelastic scattering that guarantees causality in the 
Kubo approach is in the scattering approach taken over 
by the coupling to the reservoirs. Since the electron- 
phonon relaxation is not expected to directly impede the 
overall rate of magnetic energy dissipation, we do not 
need to explicitly include it in our treatment. The en- 
ergy flow supported by the leads, thus, appears in our 
model to be carried entirely by electrons irrespective of 
whether the energy is actually carried by phonons, in case 
the electrons relax by inelastic scattering before reaching 
the leads. So wc are able to compute the magnetization 
damping, but not, e.g., how the sample heats up by it . 
According to Eq. H]), the time derivative of the energy 
reads 



E = Hcff-dM/dr = (l/7^)m G(m)m 



(2) 



in terms of the magnetization direction unit vector m = 
M/7\/s and rii = dm/dr. We now develop the scatter- 
ing theory for a ferromagnet connected to two reservoirs 
by normal metal leads as shown in Fig. [T] The total 
energy pumpin g in to both leads /^p"™^-' at low tempera- 
tures reads ll|, II4I 



jipnmp) ^ (;,/4^)Tr^^t^ 



(3) 



where S — dS/dr and S is the S-matrix at the Fermi 
energy: 



5'(m) 



r t' 
t r' 



(4) 



r and t {r' and t') are the reflection and transmission ma- 
trices spanned by the transport channels and spin states 
for an incoming wave from the left (right). The gener- 
alization to finite temperatures is possible but requires 
knowledge of the energy dependence of the S-matrix 
around the Fermi energy [12| . The S-matrix changes 
parametrically with the time-dependent variation of the 
magnetization S{t) — S'(m(r)). We obtain the Gilbert 
damping tensor in terms of the S-matrix by equating the 
energy pumping by the magnetic system ^ with the en- 



ergy loss expression ^, E ^ I 



(pump) 



E 



Consequently 
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dS dS'' 



drrii dm A 



(5) 



which is our main result. 

The remainder of our paper serves three purposes. We 
show that (i) the S-matrix formalism expanded to linear 
response is equivalent to Kubo linear response formalism, 
demonstrate that (ii) energy pumping reduces to inter- 
face spin pumping in the absence of spin relaxation in the 
scattering region, and (iii) use a simple 2-band toy model 
with spin-flip scattering to explicitly show that we can 
identify both the disorder and interface (spin-pumping) 
magnetization damping as additive contributions to the 
Gilbert damping. 

Analogous to the Fisher-Lee relation between Kubo 
conductivity and the Landaucr formula (l5[ we will now 
prove that the Gilbert damping in terms of S-matrix ([5]) 
is consistent with the conventional derivation of the mag- 
netization damping by the linear response formalism. To 
this end wc chose a generic mean-field Hamiltonian that 
depends on the magnetization direction m: 7J = H{m) 
describes the system in Fig. [TJ H can describe realistic 
band structures as computed by density-functional the- 
ory including exchange-correlation effects and spin-orbit 
coupling as well normal and spin-orbit induced scattering 
off impurities. The energy dissipation is E = (dH/dr), 
where (...) denotes the expectation value for the non- 
equilibrium state. In linear response, we expand the mag- 
netization direction m(t) around the equilibrium magne- 
tization direction moj 



m(r)= mo + u(t). 



(6) 



The Hamiltonian can be linearized as H = H^t + 
Ui{T)diH, where Hst = H{ino) is the static Hamilto- 
nian and diH = 8^11 (mo), where summation over re- 
peated indices i = x,y,z is implied. To lowest order 
E = Ui{T){diH), where 

/•OO 

{d,H) = {d,H)o + / dT'x,,{T - t')u,{t') . (7) 



(. . . )o denotes equilibrium expectation value and the re- 
tarded correlation function is 



Mr - r') = -^e{T - r') ([a,i/(r), 9,H(r')])^ (8) 



in the interaction picture for the time evolution. In order 
to arrive at the adiabatic (Gilbert) damping the magne- 
tization dynamics has to be sufficiently slow such that 
Uj{t) ~ Ujit) + {t — t)uj{t). Since m^ = 1 and hence 
rh • m = 



f 



E = id^Xiji^ -^ Q)uiU 



■J' 



(9) 



where Xijd^) ~ X_ dTXij{T)ex.\i(ioJT). Next, we use 
the scattering states as the basis for expressing the 
correlation function ([8]). The Hamiltonian consists of 
a free-electron part and a scattering potential: H = 
Hq + V^(ni). We denote the unperturbed eigenstates of 



the free-electron Hamiltonian Hq = — ?i^V^/2m at en- 
ergy e by \(ps,q{e)), where s = l,r denotes propagation 
direction and q transverse quantum number. The po- 
tential F(m) scatters the particles between these free- 
electron states. The outgoing (-I-) and incoming wave 
(-) eigenstates |V's,ij (e)) of the static Hamiltonian Hst 
fulfill the completeness conditions (V's.g (e)|^^/ i/(e')) = 
Ss,s'Sq,q'6{e — e') [lO|. These wave functions can be ex- 
pressed as \ijjs (e)) = [1 -|- Ggt Kit]|'Ps(e)), where the 
static retarded (-I-) and advanced (-) Green functions are 
G'st (^) = {e ± irj — Hst)~^ and 77 is a positive infinites- 
imal. By expanding Xiji'-^) in the basis of the outgo- 
ing wave functions lips ), the low-temperature linear re- 
sponse leads to the following energy dissipation ^ in the 
adiabatic limit 



E 



nM,u, Ui+^\d,H\4+l,) Ui:^l,\d,H\4+'> 



(10) 
with wave functions evaluated at the Fermi energy ep- 

In order to compare the linear response result, Eq. 
([TU]) . with that of the scattering theory, Eq. ([S]), we intro- 
duce the T- matrix T as ^(e; m) = 1 — 27riT(e; m), where 
f ^ V[l + G(+)f] in terms of the full Green function 
G(+)(e, m) = [e + i?7-7?(m)]~^ Although the adiabatic 
energy pumping ([5]) is valid for any magnitude of slow 
magnetization dynamics, in order to make connection to 
the linear-response formalism wc should consider small 
magnetization changes to the equilibrium values as de- 
scribed by Eq. ([6]). Wc then find 



drT^ i + VstGi-; 



(+) 



uAH 1 + G^+V, 



(11) 



into Eq. ^ and using the completeness of the scattering 
states, we recover Eq. ([T0| . 

Our S-matrix approach generalizes the theory of (non- 
local) spin pumping and enhanced Gilbert damping in 
thin ferromagncts [5|: by conservation of the total an- 
gular momentum the spin current pumped into the 
surrounding conductors implies an additional damping 
torque that enhances the bulk Gilbert damping. Spin 
pumping is an N|F interfacial effect that becomes impor- 
tant in thin ferromagnetic films [1J|. In the absence of 
spin relaxation in the scattering region, the S-matrix can 
be decomposed as ^(m) = S^{1 + a ■ m)/2 + Si{l — a- ■ 
m)/2, where <t is a vector of Pauli matrices. In this case, 
Tr(a^5') (drS)'' = Arvh^, where A = Tr[l - ReSt^j] 
and the trace is over the orbital degrees of freedom only. 
We recover the diagonal and isotropic Gilbert damping 
tensor: Gy — SijG derived earlier [5|, where 



G = jMsa 






~ ./\.r 



fl2) 



Finally, we illustrate by a model calculation that 
we can obtain magnetization damping by both spin- 
relaxation and interface spin-pumping from the S-matrix. 



We consider a thin film ferromagnet in the two-band 
Stoner model embedded in a free-electron metal 



H 



-—^^ + 5{x)V{p), 



(13) 



where the in-plane coordinate of the ferromagnet is p 
and the normal coordinate is x. The spin-dependent po- 
tential V{p) consists of the mean-field exchange interac- 
tion oriented along the magnetization direction m and 
magnetic disorder in the form of magnetic impurities Si 



V{p) =v& ■m + ^Q.a ■ SiS{p - p,), 



(14) 



which are randomly oriented and distributed in the film 
at a; = 0. Impurities in combination with spin-orbit cou- 
pling will give similar contributions as magnetic impuri- 
ties to Gilbert damping. Our derivation of the S-matrix 
closely follows Ref. [8|. The 2-component spinor wave 
function can be written as ^(x, p) = J2u Cky (a^)*&k|| (p), 
where the transverse wave function is $k|| (p) = exp(ik|| • 
P)/\/A for the cross-sectional area A. The effective one- 
dimensional equation for the longitudinal part of the 
wave function is then 

Ck||(a:) = ^fk||,k;|Ck||(0)(5(a;), (15) 



dx^ 



where the matrix elements are defined by Fkn k' — 

{2m/h^) j dp^^ {p)V{p)^\^i {p) and the longitudinal 

wave vector fcx is defined by fc^ = 2meF/fi^ — k^. For 
an incoming electron from the left, the longitudinal wave 
function is 

C ik±x<: , , , _|_ p-ik_LX t < 

Xs I ^ "kySjknS' T^ c 'k||S,k's' , J' ^ ^J 

^•""'"7^1 e'^^^ik||.,4«' ,x>0 ' 

(16) 
where s =t,i and x\ = (IjO)^ and Xl = (07l)^- Inver- 
sion symmetry dictates that t' = t and r ~ r' . Continu- 
ity of the wave function requires 1 + r = f. The energy 
pumping ^ then simplifies to /^"'"P-* — ftTr (ii'^^ /n. 
Flux continuity gives i = (1 -f iT)~^, where Fkus.k's' = 

In the absence of spin-flip scattering, the transmis- 
sion coefficient is diagonal in the transverse momentum: 
t^. = [^ — i'q^a ■ m]/(l -I- 77^), where 77^ = mv / {'h? k \^ . 
The nonlocal (spin-pumping) Gilbert damping is then 
isotropic, Gy(ni) = 5ijG' , 



9 



2z^ 
. ^(l + ,i) 



2 \2 



(17) 



It can be shown that G' is a function of the ratio be- 
tween the exchange splitting versus the Fermi wave vec- 
tor, "qp = mv /{fi?kp). G' vanishes in the limits 77^ <C 1 



orbit interaction, t 



(nonmagnetic systems) and rjp ^ I (strong ferromag- 
net). 

We include weak spin-flip scattering by expanding the 
transmission coefflcient t to second order in the spin- 

1 + toitsf - Uo«f sf j to, which 

inserted into Eq. ([5]) leads to an in general anisotropic 
Gilbert damping. Ensemble averaging over all ran- 
dom spin configurations and positions after considerable 
but straightforward algebra leads to the isotropic result 



G = G^' 



G' 



(18) 



where G' is defined in Eq. (fT7| . The "bulk" contribution 
to the damping is caused by the spin-relaxation due to 
the magnetic disorder 



^(i. 



NsS'(^'^ , 



(19) 



where Ng is the number of magnetic impurities, S is the 
impurity spin, (^ is the average strength of the magnetic 
impurity scattering, and ^ — $(^f) is a complicated ex- 
pression that vanishes when rip is cither very small or 
very large. Eq. p^ proves that Eq. ([5]) incorporates the 
"bulk" contribution to the Gilbert damping, which grows 
with the number of spin-flip scatterers, in addition to in- 
terface damping. We could have derived G*^'"*-* [Eq. (fT9|) ] 
as well by the Kubo formula for the Gilbert damping. 

The Gilbert damping has been computed before based 
on the Kubo formalism based on flrst-principles elec- 
tronic band structures [9|. However, the ab initio appeal 
is somewhat reduced by additional approximations such 
as the relaxation time approximation and the neglect of 
disorder vertex corrections. An advantage of the scatter- 
ing theory of Gilbert damping is its suitability for modern 
ab initio techniques of spin transport that do not suffer 
from these drawbacks [l6|. When extended to include 
spin-orbit coupling and magnetic disorder the Gilbert 
damping can be obtained without additional costs ac- 
cording to Eq. ([5]). Bulk and interface contributions can 
be readily separated by inspection of the sample thick- 
ness dependence of the Gilbert damping. 

Phonons are important for the understanding of damp- 
ing at elevated temperatures, which we do not explic- 
itly discuss. They can be included by a temperature- 
dependent relaxation time [3| or, in our case, structural 
disorder. A microscopic treatment of phonon excitations 
requires extension of the formalism to inelastic scatter- 
ing, which is beyond the scope of the present paper. 

In conclusion, we hope that our alternative formal- 
ism of Gilbert damping will stimulate ab initio electronic 
structure calculations as a function of material and dis- 
order. By comparison with FMR studies on thin ferro- 
magnetic fllms this should lead to a better understanding 
of dissipation in magnetic systems. 
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